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Abstract 



i c| The static three-quark potential and pseudo-potential of a pure SU(3) gauge theory are 

Oh' studied in a five-dimensional framework known as AdS/QCD. The results support the Y- 

■ ansatz for the baryonic area law. A comparison with the quark- ant iquark calculations shows 

D ' the universality of the string tension as well as the spatial string tension. We also discuss 

i-C ■ extensions to SU (N) gauge theories. 
PACS: 12.39.Pn; 12.90.+b 

(N ' 
> 
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1 Introduction 

O 

OO 

\ Heavy-quark potentials are of primary importance in studying the mass spectra of mesons and 
baryons. They have been computed in lattice simulations, and the results reveal a remarkable 
agreement with phenomenology. 1 

Until recently, the lattice formulation even struggling with limitations and systematic errors 
was the main computational tool to deal with strongly coupled gauge theories. The situation 
changed drastically with the invention of the AdS / CFT correspondence that resumed interest in 
finding a string description of strong interactions. 

In this paper we continue a series of studies [2, 3, 4] devoted to the heavy quark potentials 
and pseudo-potentials within a five-dimensional framework nowadays known as AdS/QCD. In [2], 
the model was presented for computing the quark-antiquark potential. The resulting potential 
is Coulomb-like at short distances and linear at long range. Subsequent work [5] has made it 
clear that the model should be taken seriously, particularly in the context of consistency with 
the available lattice data as well as phenomenology. In [3, 4], the models were presented for 
computing the quark-antiquark pseudo-potentials resulting from the spatial Wilson loops. The 
results obtained for the spatial string tensions are remarkably consistent with the available lattice 
data for temperatures up to 3T C . 



*Also at Landau Institute for Theoretical Physics, Moscow. 
1 For a review, see [1]. 
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The question naturally arises: What happens when these models are used for computing 
multi-quark potentials or pseudo-potentials? The multi-quark potentials have recently been the 
object of numerical studies. 2 In the case of great interest, three quark states in a SU(3) gauge 
theory, the potential is well described by a simple model [8] 

3 1 

V 3q = -a 3q — + cr 3q L min + C . (1.1) 

Kj ■> 

Here C is a constant. The quarks form a triangle with sides of lengths L^. L min is the minimal 
length of the string network which has a junction at the Fermat point of the triangle [9]. Thus, 
the potential is given by the sum of the Coulomb terms and the linear term called the Y-law. 
A remarkable fact is that <7 3q is equal to <r qq found from the quark-antiquark potential. This is 
obvious in the string picture, where the string tension is universal, but it is far from being so in 
the lattice formulation. 

The purpose of the present paper is two-fold. First, we examine the multi-quark potentials 
that may also be thought of as a further cross-check of the model [2]. Second, we make a similar 
analysis of the multi-quark pseudo-potentials. To our knowledge, there have been no studies 
(numerical or analytical) of this problem in the literature. 

The paper is organized as follows. In section 2, we discuss the multi-quark potentials. We 
begin with the three-quark potential. In this case, we demonstrate the Y-law and the universality 
of the string tension. Finally, we extend our analysis to the SU(N) case. We then go on in 
section 3 to discuss the multi-quark pseudo-potentials. Here we also demonstrate the Y-law, the 
universality of the spatial string tension, and possible generalizations to SU(N). We conclude 
in section 4 with a brief discussion of possibilities for further study. Some technical details are 
given in the appendix. 

2 Calculating the Potentials 

In this section, we will discuss static multi-quark potentials. We start with the three-quark 
potential in a SU(3) gauge theory, where equations are most elementary. In section 2.4, similar 
issues are considered in the context of a SU (N) gauge theory. 

2.1 General Formalism 

As for the quark-antiquark potential, the static three-quark potential can be determined from 
the expectation value of a Wilson loop. The baryonic loop is defined in a gauge-invariant manner 
as W 3q = ^£ a bc^a'b'c'Ui a ^3° ' with the path-ordered exponents U{ along the lines shown in 
Figure 1. 

In the limit T — > oo the expectation value of the Wilson loop is 

(WsJ-e- 2 ^, (2.1) 

where V 3q is the three-quark potential. 

2 For reviews, see [6, 7] and references therein. 
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Figure 1: A baryonic Wilson loop W^. A three-quark state is generated at t = and is annihilated at 
t = T. The quarks are spatially fixed in R 3 at points Q\, Q2 and Q3. 



In discussing baryonic Wilson loops, we adapt a formalism proposed within the AdS/CFT 
correspondence [10, 11] to AdS/QCD. 3 So, we place heavy quarks at the boundary points of the 
five-dimensional space and consider a configuration in which each of the quarks is the endpoint of 
a fundamental string, with all the strings oriented in the same way. The strings join at a baryon 
vertex in the interior as shown in Figure 2. 4 We also assume that the quarks form a triangle 



Q t U,y„0) 





Figure 2: A configuration used to calculate the expectation value of W 3 q. The quarks are set on the x-y 
plane. The baryon vertex is placed at D. Its projection on the x-y plane is Y . 

Q1Q2Q3 such that all the internal angles are smaller than 

Before proceeding to the detailed analysis, let us set the five-dimensional geometry. We 
consider the following deformation of the Euclidean AdSs [14, 2] 

sr 2 

ds 2 = R 2 w(dt 2 + dx 2 + dr 2 ) , w(r) = , (2.2) 

where dx 2 = dx 2 + dy 2 + dz 2 and 5 is a free parameter. 5 We also take a constant dilaton and 
3 For subsequent developments of the formalism, see [12, 13]. 

4 From the point of view of 10-dimensional string theory, the baryon vertex is a wrapped fivebrane whose world- 
volume is R x X, with X a 5-dimensional compact space and R a timelike curve in AdSs or its deformed version, 
[10, 11]. In AdS/QCD, it is reduced to a zerobrane (point). 

5 The value of s can be fixed, for example, from the quark-antiquark potential. 
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discard other background fields. 

The action of the system has in addition to the standard Nambu-Goto actions of the funda- 
mental strings, also a contribution arising from the baryon vertex. It is thus 

3 

S = 'y ^ Si + S vcrt , (2-3) 

i=\ 

where Si denotes the action of the string connecting the z-quark with the vertex. 

Like in the case of the quark-antiquark potential, a natural proposal for the expectation value 
of the Wilson loop is 

(VT 3q ) -e" 5 ^, (2.4) 

where 5 min is the minimal action of the system. 

Since we are interested in a static configuration, we take 

ti(n) = Ti, yi{(Ji) = mai + bi . (2.5) 
Here (rj,cjj) are worldsheet coordinates. The action of the z-string is then 

Si = Tg J dai wyja 2 + xf + rf , (2.6) 

where g = -M-,. A prime denotes a derivative with respect to <7j. 
The action for the baryon vertex is taken to be of the form 

5 vcr , = TV(r ) , (2.7) 

where V can be considered as an effective potential of the vertex. Unfortunately, the explicit 
form of V is not determined only from the 5-dimensional metric (2.2). It requires the knowledge 
of the string theory dual to QCD. We will return to this issue in the next section. 
The boundary conditions on the fields are given by 

Xi(0) = Xi , yi(0) = yi , n(0) = , ^(1) = x , y^l) = y , n(l) = r . (2.8) 

These determine the coefficients a% and b{ in (2.5). Thus, we have 



ai = y -yi, h = yi. (2.9) 

Next, we extremize the total action S with respect to the worldsheet fields Xi(ai) and r,(crj) 
describing the strings as well as with respect to x , y Q and r describing the location of the baryon 
vertex, with the following identifications: fej(l) = 5x and 5ri(l) = 5r . In doing so, we use the 
fact that there are two symmetries which simplify the further analysis. 

Since the integrand in (2.6) does not depend explicitly on <7j, we get the first integral of 
Euler-Lagrange equations 

(2.10) 



a, 2 + x? + r' 2 
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In addition, because of translational invariance along the x-direction, there is another first inte- 
gral. Combining it with (2.10) gives 

Pi = x' i . (2.11) 
Together with the boundary conditions these equations determine x% 

Xi(ai) = (x - Xi)oi + Xi . (2.12) 

Now, we extremize the action with respect to the location of the baryon vertex. After using 
(2.9) and (2.12), we get 



y Xo z^l = o, y y °~J^ = o, y - 1 + l^ (ro) = o, (2 .i 3 ) 
Ukvr+h Ukvr+ki fe./TT^ 0- 



Here 



h = (^) 2 and U = \YQi\ = y/(x - Xi f + (y - Vi f 



If we define the first integrals l{ at U{ = 1 such that Ii = w(r ) /hy/1 + and then integrate 
over [0, 1] of dai, then by virtue of (2.10) we get 



where Vi = ri/r and A = sr^. 

Now, we will compute the energy of the configuration. First, we reduce the integrals over 
Oi in Eq.(2.6) to that over rj. This is easily done by using the first integral (2.10). Since the 
integral is divergent at rj = 0, we regularize it by imposing a cutoff e. Then we replace with 
Vi as in (2.14). Finally, the regularized expression takes the form 

Its e-expansion is 

E R = ^ + 0(1). 

Subtracting the - term (quark masses) and letting e = 0, we find a finite result 



dvj 



1 + ki 



+ c, 



(2.16) 



where C stands for a normalization constant. 

In contrast to the quark-antiquark case, the potential in question is more involved. It is given 
by a set of equations. Formally, one can eliminate parameters and find E as a function of the 
/,,'s or the -Ljj's. Unfortunately, in practice it is extremely difficult. 
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2.2 A Concrete Example 

We will next describe a concrete example in which one can develop a level of understanding that 
is somewhat similar to that of the quark-antiquark case [2]. We consider the most symmetric 
configuration of the quarks, in which the triangle Q1Q2Q3 is equilateral, and specify the action 
S veit as that of a particle in a curved space. The latter implies that 



V(r )=mi?VwW, (2.17) 

where m is a parameter, which can be interpreted as a mass of a "particle". 

A consequence of the symmetry is that Y, which is the projection of D on the x-y plane, is 
nothing but the circumcenter of the triangle QiQ2Qz- The first two equations of (2.13) are now 
identically satisfied, while the last takes the form 6 



Vl + A:" 1 



K(l-A)e"2 A , (2.18) 



where k = hi and k = 



With the form (2.18), k can be explicitly computed as a function of A. Inserted in (2.14) and 
(2.16), this gives 

I = ^ p j\vv 2 e x ^(l-pv 4 e 2X ^y h (2.19) 



and 



£ = 3 0a /t 



^/>Vl_ p <;4 e 2A(l-^A *_ 1 _ v 2 



+ C, (2.20) 



A 

where p(X) = 1 - k 2 (1 - A) 2 e"\ 

The potential in question is written in parametric form given by Eqs.(2.19) and (2. 20). 7 It is 
unclear to us how to eliminate the parameter A and find E as a function of I. We can, however, 
gain some important insights into the problem from two limiting cases as well as numerical 
calculations. 

We start by noting that Eq.(2.18) makes sense only if A < 1. This means that r must obey 

r <y^. (2.21) 

So in other words, the baryon vertex is prevented from getting deeper into the r-direction. This 
gives a kind of wall which is a generic feature of confining theories. It is worth mentioning that 
the same upper bound was found in the quark-antiquark case [2] by inspecting the integral (2.19) 
at p = 1. The point is that the integral is real for A < 1. It develops a logarithmic singularity 
at A = 1 and becomes complex for larger A. A similar analysis shows that if k < 1 this holds for 
smaller p too. For k > 1 things are more subtle because of a lower bound 

r* < r . (2.22) 



6 This equation determines the location of the vertex in the r-direction, because the symmetry argument alone 
is not enough to do so. 



7 Note that these equations are reduced to those of [2] at k — and p = 1. 
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Here r* is a root of p(sr 2 ) = 0. The physical reason for this is a big mass of the "particle" that 
strengthens a gravitational force pushing the baryon vertex deeper into the interior. As a result, 
for small Z the configuration looks like a spike of height r* . 

The fact that both the equations (2.18) and (2.19) lead to the same upper bound suggests 
that neither strings nor baryon vertices are allowed to get deeper into the r-direction than l/y/s. 
This seems natural enough from the point of view of consistency. 

To complete the picture, let us present the results of numerical calculations. The parametric 
equation (2.19) predicts a characteristic form of Z, as shown in Figure 3. We see that the curves 

KM 




0.0 0.2 0.4 0.6 0.8 1.0 



Figure 3: Plots of l(X, k) for fixed k. Here s = 0.45 GeV 2 . 



behave similarly in the vicinity of A = 1 but rather differently for smaller A. This has an 
interesting effect on the form of the interquark potential, as we will see in a moment. 

Having understood the correspondence between A and Z, we can investigate the properties of 
the interquark interaction at long and short distances. 8 

At long distances the interquark interaction being independent of k is given by 

E = a 3q L min + O(l) (2.23) 

that is nothing but the desired Y-law. Here L m - m = 31 and <7 3q = egs. 

Some comments about formula (2.23) are in order. First, the form of the potential is in 
agreement with both the old flux-tube picture of hadrons and more recent calculations in lattice 
QCD [6, 7]. However, the key difference is that in our model the interaction takes place in 
the interior of five-dimensional space, whereas in those it happens on its boundary. Next, the 
tension <7 3q is the same as that found from the quark- antiquark potential [2]. Thus, the model also 
supports the universality of the string tension. Finally, to leading order the energy is the sum of 
the string contributions. A contribution from the baryon vertex appears at next-to- leading-order. 

At short distances the form of the interquark interaction depends on the value of n. We have 
the Coulomb terms 

£ = -3^ + 0(1), (2.24) 
Li 

with L a distance between the quarks, if and only if k < 1. In other cases, the energy behaves 
as E ~ const. 

8 More details are presented in the appendix. 
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Just as in the examples discussed in the AdS/CFT formulation [10, 12], in (2.24) the contri- 
bution of the baryon vertex is compatible to that of the strings. The reason for this is obvious: 
in the region of small r the metric (2.2) behaves asymptotically as Euclidean AdSs. 

As an illustration, Figure 4 shows the results of numerical calculations. We see that the value 
of re does matter at short distances, while it becomes irrelevant at large distances. 



V(GeV) 
4[ 




Figure 4: Plots of E(l, re) for fixed re. Wc set q = 0.15, s = 0.45 GcV 2 , C = 0. 

We conclude this section by making a few estimates, which might be of interest for phe- 
nomenology. 

As already noted, in AdS/CFT the baryon vertex is the wrapped fivebrane. For definiteness, 
consider type IIB string theory on AdS§ x S 5 . The D5 brane world- volume is then S" 5 x R, with R 
a timelike curve in AdS^. Now, suppose that the action of the brane is simply the Nambu-Goto 
term T5 f d 6 £,^/g, with T5 the tension. 9 In the static case this leads to the same form as (2.17). 
Moreover, we learn that re = \. 

At first sight, it seems reasonable to fix the value of re by equating a 3q and ^a: qq . 10 A little 
experimentation with Mathematica shows, however, that this results in a negative value. If we 
assume that a 3q = Ta qq , then there are positive re's for T's relatively close to 0.5, as given in 
Table 1. Note that k=\ corresponds to T = 0.27 versus T « 0.5 in lattice QCD. It is possible 

T 0.50 0.45 0.42 0.38 0.27 
re -0.11 -0.03 0.02 0.07 0.25 

Table 1: Estimates of the parameter re. 

to partially fix this discrepancy by choosing smaller re. Then the value one needs is bounded 
from above by a number of order 0.02. 

The problem is, of course, that the action for the baryon vertex we used in (2.17) is oversim- 
plified. There are extra background scalars as well as a' corrections. 11 However, there is also a 

9 The general form of the action is given by the DBI and WZ terms plus a' corrections. See, e.g., [15]. 
10 With accuracy better than a few percent, such a relation is valid in lattice QCD [8]. 

xl In contrast to AdS/CFT, where ~ -j- becomes small for large 't Hooft coupling, the phenomenological 

estimate of [2] gives %^ ~ 1. The latter suggests that a' corrections are relevant for the real world and a departure 
from the supergravity approximation is needed. 
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possibility that the string formulation being reliable at large distances provides only a qualitative 
description of the physics at short distances. 



2.3 Y-Law and String Tension 

Apparently, one of the requirements for the multi-quark model should be consistency with the 
quark-antiquark long as the fundamental strings are prevented from getting deeper into 

the r-direction, the baryon vertex, as a string endpoint, should also be prevented from doing 
so. Our next goal will be to understand the large distance behavior of the multi-quark potential 
from this point of view. 

To implement this approach, we assume that: 1) the baryon potential V(r ) is a positive and 
regular function of r such that it reaches the minimum exactly at r = \/y/s\ 2) V — > +oo as r 
tends to zero or infinity. The latter means that the vertex can't come very close to the boundary 
as well as go far away from it. 

These conditions allow us to conclude that the last equation of (2.13) makes sense only if r 
is subject to the constraint (2.22). We also learn that ki — ► as A — > 1. 

If we look at Eq.(2.14), then a short inspection shows that l{ takes large values only in the 
corner of the parameter space (ki, A) located at (0, 1). In this region ij behaves as 

l i = -\\[\H^+l-\)+0{l). (2.25) 

Note that in the symmetric case ki ~ (1 — A) 2 . Hence, the equation reduces to (A.l). 
A precisely analogous computation for E gives 

1 3 

E = --egVsYsHvki + 1 - \) + 0{1) . (2.26) 
i=i 

A contribution from the baryon vertex is of order 1, because the function V is finite at A = 1. 
Combining this with (2.25), we find the energy configuration as a function of the lengths Z, 

3 

E = a 3q J2h + 0(1) . (2.27) 
i=i 

Here o" 3q is equal to cr qq of [2], as expected from the universality of the string picture. 

However, this is not the whole story. We still need to show that Y is the Fermat point of the 
triangle formed by the quarks. The simplest way to see this is to take the limit ki — > in the 
first two equations of (2.13) which determine the location of the vertex on the x-y plane. We 
have thus 

3 1 3 1 

53 r (xo-^i)=0, £ r (yo-lfc)=0. (2.28) 

t=l i=l 

Obviously, these equations is a result of differentiating the sum Y^i=i h with respect to x and 
y . Because a solution 
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i=l i=l i=l 

determines the Fermat point, in this limit Y tends to the Fermat point of the triangle. This 
completes our derivation of the Y-law. 

2.4 The Case SU(N) 

Now we consider the case SU(N). The analysis is a little bit more complicated but the results 
are very similar. 

By analogy with section 2.1, we place N heavy quarks on the boundary of the five-dimensional 
space. Because in a generic case the quarks are not on the same plane, we need to involve the third 
spatial coordinate z. So, the quarks are fixed in R 3 at points Qi with coordinates Xj, j/j, z- L . The 
configuration of interest is constructed as the N quarks joined together by the ./V fundamental 
strings ending at the baryon vertex in the interior of the five-dimensional space whose metric 
takes the form (2.2). 

The total action of the system is given by (2.3), where the sum now runs from 1 to N . After 
choosing the gauge (2.5), the action of the i-string takes the form 

i 



Si = Tgj da, w^Jaj + xf + zf + rf . (2.30) 
The boundary conditions are given by Eq.(2.8) together with 



Zi(0) = Zi, z l (l) = z . (2.31) 

In addition to the first integrals of section 2.1, which correspond to Euler-Lagrange equations 
for the fields Xj(<7j) and rj(cjj), there is one more integral due to translational invariance in the 
z-direction. Combining it with (2.10) results in 

Pi = 4 . (2.32) 

These equations yield solutions 

Zi{<Ji) = (z - Zi)(ii + Zi (2.33) 

which coincide at the endpoints with the boundary values defined in (2.31). 

Extremizing the action with respect to the location of the vertex, we get (2.13), with the 
upper bound of summation given by N, and 

N 

that determines the location of the baryon vertex along the z-direction. Note that now Zj = 
a/ {x - Xi) 2 + (y - yi) 2 + (z - Zi) 2 . 

Like in the case iV = 3, we can compute the expressions for the lengths After doing so, we 
find that they take the forms as in (2.14) and (2.16), with the only difference: E now includes 
the contributions from the iV fundamental strings. 
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It is also straightforward to extend the analysis of section 2.3 to the case of interest. Assuming 
the same form of V, we quickly come to the similar conclusions about the bound (2.21) and large 
values of Zj. Continuing along those lines leads to the asymptotic behaviors as (2.25) and (2.26). 
Finally, we get 



As before the string tension <7 Nq is equal to that of the quark-antiquark case. 

To complete the picture, we should check that E given by the above formula is minimum. 
The essential point is to take the limit A;^ — > in the equations determining the vertex location. 
Then, from Eqs.(2.28) and 



we learn that Y is nothing but the geometric median of the set of the points Q%. This is the 
desired result which reflects the fact that E is minimum. 

3 Calculating the Pseudo-Potentials 

In this section we will investigate the temperature dependence of the spatial string tension. The 
model of interest is developed for the quark-antiquark case in [3, 4], to which the reader is referred 
for more detail. The philosophy is that the spatial string tension is determined from temperature 
dependent pseudo-potentials extracted from the spatial Wilson loops. 13 

3.1 General Formalism 

As before, we begin with the case SU(3). Let W 3q be a spatial Wilson loop in R 3 , at a fixed 
value of t, with the path-ordered exponents Ui along the lines shown in Figure 5. 



A? 




(2.35) 




(2.36) 



z 





Figure 5: A baryonic spatial Wilson loop Ws q . The points Qi are on the x-y plane. 



It is also known as the Fermat- Weber point. See, e.g., [16]. 

These are Wilson loops in hyperplanes orthogonal to the temporal direction. 
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In the limit Z — > oo the expectation value of the loop is given by 

<^ 3q )~e- z ^, (3.1) 

where V 3q is a pseudo-potential. 

Following [3, 4], we take the following ansatz for the five-dimensional geometry which turned 
out to be quite successful in the quark-antiquark case 

4 

ds 2 = R 2 w(fdt 2 + dx 2 + f^dr 2 ) , f = l- — , (3.2) 

where r T = 1/irT. Note that it is nothing but a deformation of the black hole in AdSs. 

To compute the spatial loop in question, we study a configuration similar to that shown in 
Figure 2: three fundamental strings, with the same orientation, such that the i-string begins at 
Qi and ends on the baryon vertex in the interior. The total action of the system is then given 
by the sum of the string actions and the action for the vertex. Because we are interested in the 
configuration independent of z, for the worldsheet coordinates we can choose 

Zi(Ti)=Ti, yi((Ti) = (LiCTi + bi. (3.3) 

With such a choice, the action of the i-string is then 

Si = Zsj o d<Ti w^a 2 + xf + /-iff . (3.4) 
As before, for the baryon vertex we take the action of the form 

5 vcr , = ZV(r ) . (3.5) 

The boundary conditions on the fields are given by (2.8). These determine the coefficients ctj 
and bi in (3.3). 

To find the configuration, we must extremize the total action with respect to the worldsheet 
fields Xi(ai) and ri{ai) as well as with respect to the coordinates of the vertex x , y , r . The 
underlying symmetries of the model provide a couple of first integrals that simplifies the analysis. 
Since the integrand in (3.4) does not depend explicitly on (Xj, we have an integral 

Ji = , ■ (3.6) 

yfii + x' 2 + /-irf 

Like in section 2.1, translational invariance along the x-direction yields another integral. Combin- 
ing it with (3.6) results in (2.11). Once this integral has been found, the Xj's can be determined. 
As a result, we obtain (2.12). 

Now, we extremize the action with respect to the location of the baryon vertex. At finite 
temperature, the equations (2.13) are changed as follows. There is a factor f~ l in front of k{. 
Explicitly, 

V % ;~ X \ =0, V V0 ~ y \ =0, V , 1 +lj°V(r o ) = 0, (3.7) 
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r 4 

where f Q = 1 By virtue of (3.6), the integral over [0, 1] of Oi is equal to 

' : jf - (^rt - irk** "*)"* ' (3 ' 8) 



V S(l + /o^fcj 

Here we have expressed the integration constant via the values of the fields at Cj = 1. 

Now we will present an expression for the energy (pseudo-potential) of the configuration. The 
computation is just as above. 14 At the end of the day, we have 



1 + fc% 



+c, 



(3.9) 

where C is a normalization constant. As a result, the pseudo-potential is given in parametric 
form by Eqs.(3.7)-(3.9). 

3.2 Y-law and Spatial String Tension 

Now we consider the analog of the Y-law for the pseudo-potentials. For this, we continue to as- 
sume that V(r ) is a positive and regular function with a minimum at r = \j^fs and singularities 
at the endpoints. 

For this choice of V, we can now determine the allowed values of r . A simple analysis shows 
that the equations (3.7) are consistent if r < min (l/y^, r T ) or r > max (l/y^s, r T tyl + . 
Since the latter yields complex values of Zj, it must be omitted. We have thus 

r < min(l/ v / s,r- T ) . (3.10) 

What we have found is nothing but the bound of [3, 4]. In addition to (2.22), it now states 
that neither strings nor baryon vertices may go behind the horizon (r = r T ). From the point of 
view [3], the bound gives rise to the two walls: 1) If 1/y/s < r T , then the first wall, r = 
dominates. This is interpreted as the low temperature phase; 2) If r T < 1/y/s, then the second 
wall, r = r T , dominates. This is the high temperature phase. The phase transition point is at 
r T = l/\fs. In terms of a critical temperature and the parameter s, it is written as 

T c = -V~s. (3.11) 

Having determined the allowed region for the parameter r , we are now ready to study Zj 
as a function of two variables. To this end, we look for level curves in the allowed region. 15 A 
summary of our numerical results is shown in Figure 6. The first observation is that li takes large 
values only in the vicinity of the lower right corner which lies on the A axis at A = 1 if T < T c 
and at A = 5r% if T > T c . Another important point is that at nonzero hi the length li vanishes 



14 Indeed, our previous discussion in section 2.1 corresponds to the case / = 1. 

15 As an illustration, we take ki < 2. We do not elaborate on upper bounds for the fci's, because it is irrelevant 



for our purposes. 
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Figure 6: Level curves for ij = r(A, fc,): in the low temperature phase (left) at T = T c /V2 and in the 
high temperature phase (right) at T = V~2T C . Larger values are shown lighter. Here s = 0.45 GeV 2 . 



at / = 0. This fact implies that in the high temperature phase we should keep ki/f — > to get 
large lengths. 

According to the numerical analysis, U may become large only if ki is small and r saturates 
the bound (3.10). It follows from (3.8) and (3.9), that in this limit the length and the energy 
behave as 



U = - J\/^ m (/V^ + ^/o _1 (l " fo) + /o(l - A)) P + Oil) , (3.12) 



E = -^0e A 7|^ ln(/?v^ + kif-Hl - f ) + /„(1 - A)) ? + 0(1) , (3.13) 

i=l 

where (3 2 = 4 - |fci/ _1 + - ^-f + |A(17/ - 8) - 2A 2 / . From these formulas, it is evident 
that at large distances the pseudo-potential is given by 

3 

E = a. iq J2 l * + 0(l), (3.14) 

i=l 

with the spatial tension 

if T <T C 

a 3q (J)%xp{(f) z -l} if T>T C . 



~ { q /T\ 2 (/Tc\ 2 <1 -r m ^ m (3.15) 



Here, o" 3q is the physical tension at zero temperature. The temperature dependence of cr 3q is the 
same as that of [3, 4] found in the quark-antiquark case. This suggests that the spatial string 
tension is also universal. 

Finally, taking the limit ki — > in the first two equations of (3.7), we learn that Y is the 
Fermat point of the triangle QiQ,2Q3- Thus, we have derived, in the case of pseudo-potential, 
the Y-law analogous to that of the quark potentials. 



3.3 Further Remarks 

We have here considered the five-dimensional framework for studying the spatial string tension. 
However, if one can describe the string theory construction, this gives the baryon vertex as a 
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wrapped fivebrane. In this case, the brane world-volume is R x X, with R a spacelike curve in 
the deformed AdSs. 

It is worth noting that the temperature dependence of the spatial tension is in good agreement 
with the lattice data for T < 3T C [3, 4]. Note that at higher temperatures it is determined by the 
/^-function of a pure SU(3) gauge theory [17]. Clearly, our model cannot reproduce logarithms 
associated with the running coupling. Instead, it provides a complementary description in the 
strong coupling regime. 

The above analysis can be easily generalized to the case SU (N). Compared to the discussion 
of section 2.4, the points Qi are now on the x-y plane. Hence, we needn't to involve any additional 
coordinate. The formal modification to be made is to extend the upper bound of summation to 
N. Thus, the temperature dependence of the spatial string tension is given by (3.15) that is the 
desired result [4]. 

4 Concluding Comments 

We have so far discussed the configurations including a single baryon vertex. In general, it is 
possible to consider configurations with an arbitrary number of vertices: string networks. 

To give an idea of how this works, let us consider, though only schematically, the tetraquark 
case that is the four-quark potential V 4q . We place quarks and anti-quarks at boundary points 
Qi and Qi, respectively. These points are endpoints of fundamental strings. There are two 
basic configurations, as sketched in Figure 7. The connected configuration includes five strings 
beginning at the boundary points and ending on the baryon and anti-baryon vertices in the 



interior of the five-dimensional space. The disconnected configuration, a "two-meson" state, 
includes two strings stretched between the quarks and anti-quarks. 

For the tetraquark potential, there is a well-known "flip-flop" [18]. When the quarks and anti- 
quarks are well separated, the potential is given by Coulomb terms plus a term which indicates 
the formation of the two-Y flux. On the other hand, when the quarks and anti-quarks are close, 
the potential becomes a sum of two quark-antiquark potentials. Clearly, these two possibilities 
correspond to our connected and disconnected configurations. 

From the point of view of 10-dimensional string theory, it is very easy to understand why 
the disconnected configuration dominates at small separations. The connected configuration 
includes the baryon and anti-baryon vertices that is nothing but a brane-antibrane system in 



y 




Figure 7: Configurations used to calculate the tetraquark potential. 
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10 dimensions. There is a critical separation for such a system. For smaller separations, the 
open string spectrum consists of a tachyon mode. The instability associated with the tachyon 
represents a flow toward annihilation of the brane-antibrane system. 
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Appendix 

The purpose of this appendix is to consider the two limiting cases of the equations (2.19) and 
(2.20): short and long distances. 

We begin with the case of long distances. As in Figure 3, large Us correspond to values of A 
near 1. For A — ► 1, the integrals in (2.19) and (2.20) can be evaluated to give 16 

* = ~J^ln(l-A) + 0(l), £? = -|e fl v^ln(l-A) + 0(l). (A.l) 
Now, the parameter A is easily eliminated, and one gets 



E = a 3q L min + O(l) , (A.2) 

where L min = 31 and <7 3q = egs. 

In a similar spirit, we can explore the short distance behavior of E. For k < 1, this is simple 
because we need to take the limit A — > (see Fig. 3). Expanding the right-hand sides of (2.19) 
and (2.20) in powers of A, at leading order we obtain 

1 = \h(l + O(\)), £ = -A-5(£ + O(A5)), 



where l = yj 2 F X [±, |;|;1 - k 2 } and E = 3Q y /5( 2 F 1 \\ |;1 - k 2 ] - k). Combining 
these formulas, we find 

E = -3^ + O(l) , (A.3) 
Li 

where a 3q = IqEq/^/?,. 

The case k = 1 can be treated similarly. Our previous formulas for I and E become 

Z = Z 1 A + 0(A 2 ), E = C + E^ +0(A§) , 
where l\ = l/v3s and E\ = 60^. This implies that 

E = C + ~fL^ +0(L), (A.4) 

where 7 = E\j \/37f . 



16 The integrals are dominated by v ~ 1, where they take the form J dv/^/a + b(l — v) + c(l — v) 2 . The remain- 
ing integral may be found in [19]. 
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For larger values of k, the analysis is a bit more involved because of the lower bound (2.22). 
In this case we need to take the limit A — > A* (see Fig. 3). Note that p(A*) = 0. A simple but 
somewhat tedious calculation shows that in the neighbor of A = A* the length I and the energy 
behave as 

I = (A - A*)t (I* + 0(A - A,)) , E = E* + _3J(A - A,) + 0((A - A,) 2 ) , 

with 

<° = ^(yi^^- 2 ) ■ B J=^ vA -y5(^ ert(v ^ ) - 2e " A ')- 

Here p' denotes the derivative of p at A = A*. This implies that at short distances the potential 
is given by 

E = E* + 7 L 2 + 0(L 3 ) , (A.5) 

where 7 = £7f /3(Zg) 2 . 
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